GADZOOKS! Antineutrino Spectroscopy with Large Water Cerenkov Detectors 
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We propose modifying large water Cerenkov detectors by the addition of 0.2% gadolinium trichlo- 
ride, which is highly soluble, newly inexpensive, and transparent in solution. Since Gd has an 
enormous cross section for radiative neutron capture, with ^7 — 8 MeV, this would make neu- 
trons visible for the first time in such detectors, allowing antineutrino tagging by the coincidence 
detection reaction + p e'^ + n (similarly for iJ^). Taking Super-Kamiokande as a working 
example, dramatic consequences for reactor neutrino measurements, first observation of the diffuse 
supernova neutrino background. Galactic supernova detection, and other topics are discussed. 



PACS numbers: 95.55.Vj, 95.85.Ry, 14.60.Pq 



The Super-Kamiokande (SK) water Cerenkov detec- 
tor P] has very successfully observed solar 0, atmo- 
spheric 0, and accelerator Q neutrinos. Even more sub- 
tle signals could be studied if antineutrinos were tagged 
by the detection of neutrons following + p e+ -f n 
(similarly for z?^), as the coincidence detection greatly 
reduces backgrounds from radioactivities and other neu- 
trino reactions. We propose that this could be accom- 
plished by the addition of 0.2% gadolinium trichloride 
(GdCls) to the water, and refer to the resulting detec- 
tor as GADZOOKS! {Gadolinium Antineutrino Detec- 
tor Zealously Outperforming Old Kamiokande, Super!). 
While this ability exists in < 1 kton detectors, it has 
never been achieved on the scale of the 22.5 kton fiducial 
mass of SK, nor in a light water detector. The proposed 
use of GdCla dissolved in water is also novel. We discuss 
SK as a concrete example to show feasibility, but the 
proposed technique has much more general applicability. 

Neutron Capture on Gadolinium. — Neutrons in 
water quickly lose energy by collisions with free protons 
(and oxygen nuclei); once thermal energies are reached, 
the neutron continues to undergo collisions, changing its 
direction, but on average not its energy, until it is cap- 
tured j5| . The cross section for thermal neutron capture 
on natural Gd is 49000 barns, compared to 0.3 barns on 
free protons 0. With the proposed 0.2% admixture by 
mass of GdCls in water, 90% of neutron captures are 
on Gd (8 MeV gamma cascade), 0.2% on CI (8.6 MeV 
gamma cascade), and the rest on H (2.2 MeV gamma, 
not detectable in SK). After thermalization, capture oc- 
curs in about 20 ^s (about 10 times faster than in pure 
water) and about 4 cm; both arc slightly increased by 
the pre-thermalization phase. Compared to typical time 
and distance separations between events in SK, as well 
as the position resolution, these are exceedingly small. 

Since SK was first proposed, the price of the water 
soluble salt GdCls has fallen three orders of magnitude, 
most of that in the past few years. The current price of 
99.99% pure GdClg • XH2O from the Stanford Materials 
Corporation is about $3/kg. The properties of GdCls are 
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discussed in more detail below. Neutron capture on Gd 
leads to 8 MeV shared among 3 — 4 gammas. In scintilla- 
tors, where Gd is used routinely, the summed gamma en- 
ergy is relevant. However, in water what matters are the 
electrons Compton-scattered above the Cerenkov thresh- 
old by relatively hard gammas. The detectable light fol- 
lowing neutron capture on Gd (in thin foils) in possible 
discrete counters was carefully simulated 01 for the SNO 
heavy water Cerenkov detector. The equivalent single 
electron energy was found to peak at about 5 MeV, and 
range over 3 — 8 MeV. This spread reflects the intrinsic 
variation in the gamma cascades and the detector energy 
resolution (the simulation had 5 photoelectrons per MeV 
of electron energy, compared to 6 in SK-I). 

At energies as low as 3 MeV, background singles rates 
are much too high to detect neutron captures in isolation, 
as required in SNO. However, using the detection reac- 
tion j/g -f p — > n and requiring a very tight time and 
position coincidence between the positron and the follow- 
ing neutron will allow detection of nearly all captures on 
Gd. Due to continuing improvements in online data ac- 
quisition and filtering, it is expected that SK-HI, which 
will begin operations in mid-2006, will trigger at 100% ef- 
ficiency at 3 MeV and above, with good trigger efficiency 
down to 2.5 MeV. The rate of accidental coincidences is 
vanishing, since the number of candidate solar neutrino 
events is ~ 1/yr/ton in SK. Gamma cascades also pro- 
duce more isotropic light than other backgrounds, which 
aids in their identification. 

Neutron Rates in SK. — The ambient neutron rate 
in the SK fiducial volume is unknown, as neutrons and 
their captures on free protons are invisible. Though a 
high background neutron rate could be tolerated for the 
coincidence signal -\- p ^ e"*" -I- n, we do not want to 
compromise the 5 — 20 MeV solar neutrino singles signal 
by making neutron captures visible. There are ~ 10 sig- 
nal events per day along the solar direction, and ^ 100 
background events per day over all directions. After ex- 
amining a wide variety of potential sources of ambient 
neutrons, we find that they are not a problem. 
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Cosmic ray muon spallation in SK produces a very 
steeply falling spectrum of ~ 10^ neutrons/day; their 
production is of interest Spallation also produces 
beta-unstable daughters that form a significant back- 
ground for solar neutrinos; after cuts following muons, 
< 100 spallation betas/day survive. Since these cuts veto 
events within seconds and meters of muon tracks, they 
will be much more effective on neutrons, which capture 
in much less time and distance. Neutrons from spalla- 
tion events in the surrounding rock are stopped by 4.5 
m of water shielding surrounding the fiducial volume. 
Our conclusions are supported by data from KaniLAND, 
which is adjacent to SK and has similar shielding 9]. Dis- 
solved Gd would also be present in the optically isolated 
outer detector, and the neutron capture rate there could 
be as high as 10'^ Hz, due to neutrons from the surround- 
ing rock 0. However, captures on Gd would not produce 
enough light to trigger the sparse outer detector. 

The rate of i^e + P — > e"*" + interactions in SK from 
nuclear reactors is ~ 30/day. For the majority of events 
the positron would be detectable in coincidence with the 
neutron (even without that, the excess singles rate due 
to neutron captures could be detected). For the lowest 
energy reactor antineutrinos, as well as those originating 
from U/Th decays in Earth, the positron energy is too 
low to trigger SK, and the neutron would appear in iso- 
lation. The rate due to U/Th decays is expected to be 
~ 4/day. Atmospheric neutrino neutral-current events 
may be otherwise invisible if only a neutron is scattered, 
a relatively common event, at a rate of ^ 2/day. 

At the present U/Th/Rn concentrations in SK, < 1 
neutron/day in total is produced by the following pro- 
cesses: spontaneous fission of ^^^U; (a, n) reactions on 
^H, ^^O, and ^*0; and ^H(7,n); estimated by scaling 
from SNO results 0|. The GdCla additive must meet 
radiopurity standards about 10^ times less stringent than 
for the SK water. Initial test samples, for which no spe- 
cial care was taken, were measured by mass spectros copy 
to have p^^U] ~ IQ-** g/g and p^Th] ~ 10'^° g/g 
Assuming secular equilibrium in the decay chains, the 
beta and neutron rates in SK would remain similar to 
present values if the samples were purified by a factor of 
100. The Palo Verde experiment obtained Gd 10 times 
more radiopure than our initial samples, and the SK 
water system reduces U/Th/Rn by orders of magnitude 
from the original mine water. We are confident that the 
desired radiopurity is easily obtainable (as was Ref. 01). 

Natural Gd contains 0.2% ^^^Gd, which alpha decays 
(Ti/a = 1014 yr, T„ = 2.1 MeV) @. With 100 tons 
of GdCla in SK, the decay rate is ~ lO^o/day These 
alphas are invisible in SK, but their introduction may 
initiate i^O(a,n) and ^^0{a,n) reactions. Using the al- 
pha stopping power and the measured cross sections > 
the neutron production rate is ~ 1/ day. Lanthanide con- 
taminants (< 10""*) and their decays can also be ignored. 

Reactor Neutrinos. — The total + P ^ e'^ + n 



rate in SK from reactors can be scaled from the Kam- 
LAND rate |^, and is ~ 30/day after oscillations. The 
threshold for solar neutrinos in SK-I was — 5 MeV, 
and SK-III should be even better due to a lowered trigger 
threshold and much-improved offline reconstruction algo- 
rithms currently being evaluated in SK-II. This was the 
analysis threshold for single events. The trigger is effi- 
cient to much lower energies (as low as 3.5 MeV in SK-I), 
where radioactivity backgrounds overwhelm the solar sin- 
gles rate. However, for the coincidence signal, it should 
be possible to identify real events as low as Eg = 2.5 
MeV, covering most of the reactor spectrum. The present 
KamLAND analysis threshold is Eyis — + 2me = 2.6 
MeV, corresponding to Ee = 2.1 MeV in SK. Measure- 
ment of the positron energy determines the neutrino en- 
ergy, since Ei, ~ Ee + 1-3 MeV, and additionally there is 
a weak directional correlation |5j . 

GADZOOKS! would have the advantage of much 
larger statistics, with about 50 times more fiducial mass. 
The very high rate would allow the fiux to be moni- 
tored on yearly basis with about 1% statistical error, 
likely allowing new tests of neutrino oscillation param- 
eters as reactors at different distances go through on/off 
cycles. The expected reactor spectrum is shown in Fig. 1. 
The energy resolution in SK is about 6 times worse than 
in KamLAND, so that spectral distortions (not shown) 
will be smeared, though high statistics may still reveal 
them. Resolution is why the positron spectrum extends 
to Ee = 12 MeV, even though the neutrino spectrum 
plummets beyond E^ ~ % MeV. 

In KamLAND, fast neutrons from the walls can scat- 
ter protons (producing scintillation light) and then cap- 
ture, mimicking the coincidence signal; in SK, the struck 
protons are invisible. Muon spallation produces three 
isotopes (*He, ®Li, and ^^Li) that beta decay to excited 
daughter states that decay by neutron emission 0. In 
KamLAND, spallation cuts control these to be ^ 2% of 
the reactor signal. Spontaneous fission of ^ssy and at- 
mospheric neutrino neutral-current interactions can both 
produce multiple neutron events that can mimic the re- 
actor signal, but the rates are small. 

Diffuse Supernova Neutrino Background. — 
From a supernova in the nearest galaxy, Andromeda, SK 
would detect ^ 1 event. While more distant galaxies 
have lower fluxes, there are many more of them, and 
the flux from all previous core-collapse supernovae in the 
universe may be detectable. We refer to this as the Dif- 
fuse Supernova Neutrino Background (DSNB), and not 
the more common "relic supernova neutrinos," since the 
latter causes confusion with "relic" big bang neutrinos. 
All neutrino flavors are produced, but Ue is the easiest 
to detect, via Ue + p ^ e'^ + n on free protons. The 
strongest limits are from SK, based on the electron spec- 
trum above 18 MeV (at 18 MeV, the measured rate is 
- 1/22.5 kton/yr/MeV) |li|. The DSNB spectrum is 
a convolution of a single supernova spectrum with the 
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supernova rate as a function of z, with neutrino energies 
redshifted as E^/{l + z). The Kaphnghat, Steigman, and 
Walker (KSW) model pushed uncertainties in the direc- 
tion of producing the largest reasonable DSNB flux |0| . 
In the relevant energy range, E = 10 — 30 MeV, other 
models have nearly the same shape but differ mostly in 
normalization; also, several uncertainties are minimized. 
The supernova rate is reasonably known in the relevant 
range z S. 1 (where it rises by ^ 10 over the z = 
rate) [la, llal- Uncertainties on cosmological and neu- 
trino oscillation parameters no longer play a significant 
role 01, the latter especially if realistic neutrino tem- 
peratures are used. The DSNB detection cross sec- 
tion may be treated at lowest order at present. 

In Fig. 1, we show a range of DSNB spectra. The up- 
per edge of the band is set by the SK limit T^l (0.6 of 
KSW), and the lower edge by modern mode ls llQl (0.2 of 
KSW). The background-limited SK search [lj| will gain 
the required factor of 3 in sensitivity in about 40 years. 
In GADZOOKS!, this sensitivity would be available im- 
mediately. Requiring neutron detection would dramat- 
ically lower the backgrounds below 18 MeV, where the 
spallation beta singles rate rises rapidly. As the thresh- 
old is lowered, the atmospheric neutrino backgrounds fall 
and the DSNB signal rises. We calculate that the atmo- 
spheric [Ts] backgrounds can be reduced by ~ 5 from the 
measured rates |l4j by rejecting events with a preced- 
ing nuclear gamma [ly| or without a following neutron. 
Further rejection (not shown) is likely possible by re- 
quiring a small position separation between prompt and 
delayed events, since DSNB events produce much less en- 
ergetic neutrons. The number of DSNB events expected 
in GADZOOKS! is about 2 - 6 per year above 10 MeV. 
Uncertainties smaller than the corresponding Poisson un- 
certainty can be ignored. 

Detection of the DSNB would be an extremely impor- 
tant scientific milestone. It could be the first detection of 
neutrinos from significant redshifts z "^1, and the second 
detection of supernova neutrinos. With the exception of 
SN 1987A in the Large Magellanic Cloud, a close com- 
panion of our Galaxy, neutrinos have never been detected 
from farther than the Sun. The DSNB flux is propor- 
tional to the rate of all core-collapse supernovae, includ- 
ing optically dark "failed" supernovae that collapse to 
black holes. The DSNB spectrum shape would also pro- 
vide a crucial calibration for numerical supernova models, 
since in relatively few years, the sample of neutrinos from 
SN 1987A could be surpassed 

Galactic Supernova. — Supernovae in our Galaxy 
are expected about 3 times per century, and SK would 
observe ~ 10^ events at a typical distance of 10 kpc. The 
ability to cleanly identify the dominant D,, -\- p ^ ~\- n 
events would be extremely important for studying the re- 
maining reactions, notably O — > F, which 
is exquisitely sensitive to the temperature and hence 
neutrino mixing |2o| . Hundreds of Ve events could be ob- 
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FIG. 1: Spectra of low-energy Ve + p ^ + n coincidence 
events and the sub-Cerenkov muon background. We assume 
full efficiencies, and include energy resolution and neutrino os- 
cillations. Singles rates (not shown) are efficiently suppressed. 



served, far more than in any other detector. The neutral- 
current events on ^^O that lead to gamma and/or neutron 
emission would be much better identified, of key impor- 
tance for measuring the Vf^/vr temperature |2]| . Even in 
the forward angular cone, inverse beta events dominate 
neutrino-electron scattering events ^2^. Isolating those 
events would be very useful: it would help detect the neu- 
tronization burst v^. events; it would improve the pointing 
to the supernova by a factor of about 2, down to about 
2° 22]; and it would allow better spectral studies. Using 
timing information alone, SK could immediately recog- 
nize a supernova as genuine by the unique time structure 
of the events: almost all events in pairs separated by 
tens of microseconds, much shorter than the separation 
between subsequent neutrino interactions. Neutron de- 
tection would improve the ability to study bursts out to 
very late times, or to detect faint bursts. 

Other Physics. — The solar De flux is < 1% of the 
predicted flux |2^ . Requiring a neutron coincidence 
would greatly reduce backgrounds, and the sensitivity 
would be better than about 0.01%; the weak directional 
information |^ may allow even better sensitivity. For 
atmospheric and accelerator neutrinos, the ability to de- 
tect neutron captures delayed from the neutrino inter- 
action would shed new light on the hadronic flnal state. 
This would be useful for separating (especially sub-GeV) 
neutrinos and antineutrinos (which preferentially eject 
protons and neutrons, respectively) and for probing the 
type of neutrino interaction. For accelerator neutrinos, 
neutrino-neutron elastic scattering will be a signiflcant 
new neutral-current channel. For proton decay, neutron 
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detection may reduce atmospheric neutrino backgrounds, 
as well as aid the study of bound-nucleon decay modes 
(including invisible modes) in which the final nucleus de- 
cays with an emitted neutron. 

R&D on GADZOOKS! — Many of the properties 
of GdCla are already well known. It is naturally highly 
water soluble, in concentrations up to 50%. In contrast, 
an elaborate chemistry is needed to suspend even modest 
amounts of Gd in liquid scintillator. GdCla is stable, non- 
reactive, and easy to handle in bulk quantities. Years of 
human and animal studies have demonstrated its non- 
toxicity |24j . and it is commonly employed as a safe and 
effective contrasting agent in human MRI subjects. Still, 
because of the novel use we are proposing for GdCls, 
some new, specialized information is desired. 

In our calculations, we treated neutron capture on Gd, 
but otherwise ignored the effects of GdCla on the detector 
(note that Haxton discussed adding 5 — 10% admixtures 
of various other salts to SK as solar neutrino targets, but 
not for neutron capture |25|). Three key R&D issues are 
currently being addressed |26| . First, while preliminary 
measurements indicate that the light absorption length 
over the Cerenkov frequency range in Gd-loaded SK wa- 
ter remains > 100 m [llj, its complete light attenuation 
characteristics (including scattering) over these long dis- 
tances will be evaluated. Second, the physical effects 
of GdCls on detector components will be investigated 
through extended exposure of samples of these materi- 
als, as well as through the use of established accelerated 
aging techniques. Finally, an operational replica of the 
SK water filtration system will be constructed in order to 
determine the most effective method of filtering the water 
of other impurities without incurring unacceptable losses 
or concentration variations of the 0.2% GdCls solute in 
the process. 

Conclusions. — We propose GADZOOKS!, a large 
water Cerenkov detector that would allow neutron de- 
tection by radiative capture on a dissolved gadolinium 
salt. Extensive research and preliminary R&D, briefly 
reported here, strongly support the feasibility of this 
technique, and a program of focused R&D has just been 
funded M- The new ability to detect neutrons in a large 
light water detector would allow the clear identification of 
Ve by the coincidence detection reaction De+p + n 

(similarly for v^). An entirely new program of antineu- 
trino spectroscopy would be opened, with important im- 
plications for reactor, solar, supernova, atmospheric, and 
accelerator neutrinos. The prospects for the first detec- 
tion of the diffuse supernova neutrino background are 
particularly exciting. This could be the first detection of 
neutrinos from cosmological distances (and certainly the 
largest L/E), and the second detection of supernova neu- 
trinos. In just a few years, the yield from SN 1987A could 
be surpassed. Unlike previous neutron detection tech- 
niques, ours is scalable at reasonable expense. Megaton- 
scale water detectors with GdCla could observe hundreds 



of DSNB events per year, allowing stringent tests of the 
black hole formation rate and supernova neutrino spec- 
tra. Other physics topics would enjoy similar benefits. 
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